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Abstract. We present results of an analysis of all 480 ks of Rossi X-ray Timing Explorer Proportional Counter Array data 
obtained from 17 May 1998 to 11 October 1998 on the luminous low mass X-ray binary GX 13+1. We analysed the spectral 
properties in colour-colour diagrams (CDs) and hardness-intensity diagrams (HIDs) and fitted the power spectra with a multi- 
Lorentzian model. GX 13+1 traces out a curved track in the CDs on a time scale of hours, which is very reminiscent of a 
standard atoll track containing an island, and lower and upper banana branch. However, both count rate and power spectral 
properties vary along this track in a very unusual way, not seen in any other atoll or Z source. The count rate, which varied by a 
factor of ~ 1.6, along a given track first decreases and then increases, causing the motion through the HIDs to be in the opposite 
sense to that in the CD, contrary to all other Z and atoll sources. Along a CD track, the very low frequency noise uniquely 
decreases in amplitude from ~5 to ~2% (rms). The high frequency noise amplitude decreases from ~4% to less than 1% and 
its characteristic frequency decreases from ~10 to ~5 Hz. The 57-69 Hz quasi-periodic oscillation (QPO) found earlier is also 
detected, and no kHz QPOs are found. In addition the entire track shows secular motion on a time scale of about a week. The 
average count rate as well as the amplitude of the very low frequency noise correlate with this secular motion. We discuss a 
possible explanation for the peculiar properties of GX 13+1 in terms of an unusual orientation or strength of a relativistic jet. 
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1. Introduction 

Low mass X-ray binaries (LMXBs) containing a neutron star 
can be divided into two subclasses, atoll and Z sources, 
based on their correlated X-ray timing and spectral properties 
(Hasinger & van der Klislll989l) . In the colour-colour diagram 
(CD) LMXBs tend to trace out a well-defined, fully connected 
1 -dimensional track on a time scale of hours to days for Z 
sources and days to weeks for atoll sources. 

Z sources trace out a Z-shaped track in the CD, where the 
three branches of the Z from top left to bottom right are called 
the horizontal branch (HB), the normal branch (NB) and the 
flaring branch (FB). The (sharp) turns in the track separat- 
ing these branches are called vertices. The power spectral fea- 
tures found in Z sources are very low frequency noise (VLFN), 
which has an approximately power law shape and is observed 
up to a few Hz, a band limited noise (BLN) component with 
a characteristic frequency of 2 to 20 Hz and various differ- 
ent quasi periodic oscillations (QPOs) with a Lorentzian shape 
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and f requenci es that range from a few Hz to more than a kHz 
Jvan der Klisll995tlvan der KlisfeOOOl) . 

These power spectral features occur correlated to the posi- 
tion of the source in the Z track (see Fig. QPOs with fre- 
quencies of 13-55 Hz are found on the HB (HBOs), and with 
frequencies of 4-7 Hz on the NB (NBOs); in some sources the 
latter continue onto the FB where the frequency increases up 
to ^20 Hz. BLN with a fractional root-mean-square amplitude 
(rms) of up to ~8% is found along the HB and disappears as 
the source moves onto the NB. There is a strong tendency of 
the amplitude of the VLFN to increase from less then 1% to 
~10% along the whole track from the HB to the FB. The flux 
increases from left to right along the HB, stays approximately 
constant or slightly decreases along the NB, and usually in- 
creases again on the FB. The mass accretion rate (M) is thought 
to increase along the H B via the NB to the FB o n time scales 
of hours to days ( lHasinger & van derKlisl fl989). although on 
longer tim e scales the relation with M may be mo re compli- 
cated (e.g. iKuulkers et alJl994lHoman et alJ2002t) . 

Atoll sources trace out a track that generally has a C or 
U-shape, consisting of an island branch, which often takes the 
shape of one or more isolated patches because of observational 
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Fig. 1. Typical CDs (top) of Z (left) and atoll (right) sources, and power spectra (bottom) corresponding to their different branches 
in the CD. The power spectra of the HB and the island state (IS) were shifte d upwa rds by a factor of 10 4 and those of the NB and 
lower banana (LB) by a factor of 10 2 . The figure was taken from lWii nandsl J200 ll) . 



windowing (usually in the upper left region of the CD), and a 
banana branch (extending from the lower banana (LB) -lower 
left- to the upper banana (UB) -upper right). VLFN, BLN and 
QPOs simil ar to those in Z sources also occur in atoll sources 
(see Fig.Q] IWiinandsl200ll) . Along an atoll-source track from 
the island to the lower, middle and then to the upper banana, 
the rms amplitude of the VLFN increases from ~2 to 5%, the 
characteristic frequency of the BLN, which som etimes consists 
of several components Jvan Straaten et all20 02). first increases 
and then decreases within a range of ^2 to ^70 Hz, and its 
rms decreases from ^20% to l ess than 1%. The count rate in- 
creases by a factor of ~2-10 jHasinger & van der Klislll989l) 
going from the island state to the banana state (except for the 
transient sources, see below). The frequencies of the various 
QPOs increase in the same sense along the track in Z and atoll 
sources. 

GX 13+1 is a LMXB containing a neutro n star, since 
X-ray bursts have been fo und from this source jFleisch man 
Il985t iMatsuba et al.lll995l) . Together with GX 3+1, GX 9+1 
and GX 9+9, GX 13+1 forms the subclass of persistently 
bright atoll-sources (the "GX atoll so urces") which were ob- 
served to be in the banana state by Hasin ger & van der Klis 



1989). In GX 13+1 and GX 3+1, two branched structures 



gram (HID) (Stella et alll985llLewin et al.ll987HSchulz et alJ 


1989; iHoman et all 


1998; 


Muno et al.l 


2002). In luminos- 



ity, these sources are intermediate between the very lu- 
minous Z source s and the weaker remaining atoll sources 
(For d et alJl2000l) . Contrary to the other atoll sources and the 
Z sources, these GX atoll sources have so far not shown 
kHz QPOs. Although GX 13+1 has been classified as an 
atoll source based on u pper-banana branch phenomenology 
( Hasinge r & van der Klisl[l989l) . among the atoll sources it is 
the source that shows X-ray properti es closest to those of the 
Z sources (VLFN amplitude: Hasin ger & van der Klisl ll989. 
CD track: Muno et al. 2002, ~65 Hz QPOs: Homan et al. 
1998). The radio luminosity of GX 13+1 is similar to that 
of the Z sources, whereas the other atoll s ources are weaker 
JGarcia et alJll988l iFender & Hendrvll2000l) . From infra-red 
specfioscopv lBTndvopadlivav*eTaLr i 19991) conclude GX 13+1 
might be a Z source, based on the nature of its companion sug- 
gesting a long orbital period. 

High resolutio n spectrosc opic observations o f GX 13+1 
with ASCA dlJeda etall200ll) and XMM-Newton JSidoli et alJ 
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2002J, have revealed the presence of narrow absorption lines, 
similar to those in the strong galactic jet sources GRS 
1915+105 and GRO J1655-40 and the dippers MXB 1659-29 
and X 1624-49. 

For some Z sources it is known that the position of the 
track they trace out in the CD drifts on a time scale of weeks, 
with shif ts in the co lours of up to ~10% (iKuulkers et alJl 994: 
iKuulkers et aljfl99^) : this effect is called secular motion. It is 
also seen in the atoll source 4U 1636-53 JPrins & van der Klisl 
1 19971 iDi Salvo et"ai1l2003l) . but in th e CD it is not as s trong 
as in some Z-sources. In GX 13 +1. iMu no et al.l J2002I) also 
detected secular motion (see also Ho man etalJll998h . These 
shifts are usually larger in the HIDs compared to the CDs due 
to the shift of the track being strongly correlated with X-ray 
luminosity (L x ). The power spectral properties of all LMXBs 
seem to depend primarily on the position (S, Hasin ger et alJ 
Il990l) of the source in the atoll or Z track and not on L x . L x 
shows a strong correlation (usually positive, sometimes nega- 
tive) with this S on short time scales (hours to days), but this 
correlation can be completely absent on long time scales. It has 
long been assumed that this is because S is a measure of the 
accretion rate M, with L x not proportional to M due to uncer- 
tain causes such as, e.g., beaming, but more recently a scenario 
has been suggested in which the inner radius of the accretion 
disc (ri nncr ) determines the power spectral properties, while at 
the same time L x is deter mined by the total mass accretion plus 
perhaps nuclear burning ( van derKlisl200l1) . The decorrelation 
of L x from the power spectral properties occurs in this scenario 
through energy release that does not originate from accretion 
through the inner disc, but does respond to inner disc count 
rate changes in a time-averaged sense. 

Recentlv. lMuno et al.ld2002l) and lGierlihski & Dond(l2002l) 
suggested that the distinction between atoll and Z sources 
might be an artifact of incomplete sampling, because the atoll 
sources that cover a wide enough range in count rate (the tran- 
sient sources Aql X-l, 4U 1608-52 and 4U 170 5-44, each with 
a flux range of more than a factor of 10, see also lBarret & O live 
120021) also exhibit roughly Z-shaped tracks with at the low- 
est luminosities a 'horizontal branch'. However, the differ- 
ences in flux ranges and time scales on which these tracks 
are tr aced out in the CDs (as also noted by IMuno etall 
120021) as well as their to pology and power spectral behaviour 
( van Stra aten et alJ2003l) . argue against any simple unification 
of these two classes. Being a persistently bright source with a 
flux range of a factor of ~2, GX 13+1 is not among the sources 
that has been observed to show such an atoll source 'horizontal 
branch' . 

In this paper we discuss the X-ray colour, flux and tim- 
ing properties of GX 13+1 and their mutual dependence. We 
find highly unusual power spectral - colour-colour diagram be- 
haviour, which does not fit in with either Z or atoll behaviour 
as well as very unusual secular motion. In Sect. l4.4l we suggest 
a possible explanation for the unusual properties we find in GX 
13+1. 



2. Observations 

We used all publicly av ailable Rossi X-ray Timing Explorer 
(RXTE. iBradt et aljfl993l) Proportional Counter Array (PCA, 
LFahodaetalJll99d) data of GX 13+1 in the third RXTE gain 
epoch, a total of 44 observations (listed in Table Q ranging 
from May 17 to October 10, 1998 and lasting about 480 ksec 
in total. Individual observations typically last ~10 ksec, in- 
terrupted several times by South Atlantic Anomaly passages 
and Earth occultations. The uninterrupted parts of observations 
are normally ~3 ksec in duration. In 33 observations the five 
Proportional Counter Units (PCUs) 0-4 were all active, but in 
1 1 observations PCU 3 and/or 4 were inactive for part or for the 
entire observation. The 2-60 keV count rate for 5 PCUs varied 
between ^3300-5500 cnts s _1 . For all observations Standard 
2 data are available which have a time resolution of 16 s in 129 
energy channels (covering the 2-60 keV PCA range), as well as 
high time resolution Single Bit and Event data with a time res- 
olution of 1/4096 s or better, also covering the 2-60 keV PCA 
range. 

Standard 2 data of PCUs 0, 1 and 2 (which were always ac- 
tive) were used to create CDs and HIDs with 16 s data points. 
The data were background corrected using the Bright Source 
Model for the RXTE PCA 1 . The soft colour (SC) was defined 
as the count rate ratio between PCA Standard 2 spectral chan- 
nels 7-14 and 3-6 (3.6-6.5 and 2.2-3.6 keV respectively) and 
the hard colour (HC) as the ratio between channels 21-40 and 
15-20 (8.7-16.0 and 6.5-8.7 keV). The intensity was taken to 
be the 2-60 keV (channel 1-129) count rate. 

We checked the colours for instrumental drifts caused by 
changes in the response of the RXTE PCA by analysing 4 Crab 
observations covering the entire observation period of the data 
we used. We created CDs of the Crab observations using the 
same bands as we used for the GX 13+1 data and calculated 
the average colours. The shifts of Crab in the CD are below 
1% in SC and below 0.5% in HC. Dead time effects result in 
shifts of the data points in the CD relative to one another of less 
than ~0. 1 % in both HC and SC. Both the effect of instrumental 
drifts and of dead time are negligible compared to the effects 
reported below. 

3. Analysis and results 

3. 1. The colour-colour and hardness-intensity 
diagrams 

On cursory examination, the CD of our data appears to exhibit 
the lower part of an atoll track or perhaps the NB/FB part of a Z 
track, traced out on a time scale of hours to days (Fig. [2} with a 
rather sharp vertex. This track shifts on longer time scales, re- 
sulting in a CD with overlapping tracks. Similar shifts are also 
observed in the HIDs. We divided the data into 15 different 
tracks by adding observations together in chronological order 
for at most one week, unless the track shifted before then (see 
Table [Q. A shift in this context was defined as motion of the 
track over a colour-colour interval more or less perpendicular 
to the track, significant compared to the track width. Clearly, 

1 http://rxte.gsfc.nasa.gov/docs/xte/whatsnew/calibration.html 
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Fig. 2. The colour-colour diagrams of tracks 1 to 12 (black points); the tracks that show a vertex are plotted in the left plot and 
those that do not in the right. For comparison the colour-colour diagram of track 6 is plotted (grey points) in every frame. Notice 
the shifts of the track, mainly in the hard colour. For each point we used a 16 s average spectrum to calculate the SC (3.6-6.5 
keV/2.2-3.6 keV) and HC (8.7-16.0 keV/6.5-8.7 keV). Typical errors are 0.04 (top left of the CD) up to 0.06 (top right of the 
CD) in the SC and around 0.013 in the HC. 



a shift parallel to the tracks can be overlooked in this way. For 
that reason data from our 15 tracks could still be to some extent 
mixed up. However, it is not likely that this happened very of- 
ten. The reason for this is, that, as shown below, the shifts tend 
to occur along one particular diagonal in the CD. The tracks do 
not usually run parallel to that diagonal. 

Tracks 3, 10a,b,c and lla,b were excluded from further 
analysis, because they do not show a vertex. Because of this 
lack of a vertex we cannot with confidence compare the posi- 
tions of these tracks in the CD with those of the others. 

For each of the remaining 9 tracks we defined a parameter 
S a that indica tes the position along the track, in a similar way 
as e.g. Dieters & van der Klis (2000) for the Z source Sco X- 1 . 
This was done by approximating the track with a spline onto 
which the data points are projected along a vector defined by 
their errors in HC and SC (see Fig.|5Jl- Different from e.g Sco 
X-l, GX 13+1 only shows one vertex in the CDs instead of two. 
We chose the value of 5 a at this vertex to be equal to 2, and 
fixed the arc length along the track by setting S a =l at a point 
on the top branch at a fixed HC and SC interval (AHC=0.3, 
ASC=0.2) above and to the right of the vertex. The locations 
of the vertices (the middle of the sharp turn) were estimated by 
eye. To quantify the shift of each track in the CD we defined 
a parameter S s hm as the distance of the track's vertex to the 
po int HC=0.46 and S C=3.07 (the vertex of track 1): S^ift = 
^(AHC 2 + ASC 2 ). 



In Fig. |2] we show the CDs of all 15 tracks in black su- 
perimposed on track 6 (in grey) for reference. As mentioned 
above GX 13+1 shows a pattern very reminiscent of a stan- 
dard atoll track, with, however, a rather sharp vertex. The 
tracks are shifted relative to each other; the vertex locations 
are presented in Fig. 0] and Tabled The shifts are consistent 
(X 2 /dof = 8.3/7) with being on a diagonal line in the CD. 
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Fig. 3. The CD of track 6 and the spline used to calculate S a . 
The locations on the curve where S & = 1 and S a — 2 are 
indicated; the arrow indicates the direction of increasing S a . 
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Fig. 4. The vertex locations of the tracks with their track num- 
ber and the best linear fit to the points. The vertices are con- 
sistent with being on a line of increasing hard and soft colour 
(X 2 /dof = 8.3/7), spanning a 7% increase in SC and a 43% 
increase in HC. 



Representative HIDs of track 1 are shown in Fig.|5] together 
with the CD. The colour-colour points have different symbols 
based on their location in the CD to simplify tracking the differ- 
ences in the shape of the track between the diagrams. Contrary 
to what is seen in other sources the track is flipped in the HIDs 
compared to the CD. 

3.2. Dependence of the count rate on S a and S shiit 

In Fig. [6] we show the 2-60 keV count rate plotted ver- 
sus S a for each track. Contrary to what is normally seen 
in atoll sources on these timesca les llHasinger & van der Klisl 
Il989tbi Salvo et alJl200lLl2003l) . the count rate depends non- 
monotonically on S a , with usually a minimum somewhere in 
the range S a =2-3.4. Such a minimum is seen at th e NB-FB ver- 
tex (S a =2) in some Z sources (see, e.g. iDieters & van der Klisl 
200(i l.Tonker et alJ200oT:lHoman et alJ2002h . but there the min- 
imum always occurs at the vertex, while in GX 13+1 its po- 
sition drifts between S a values of 2 and 3 depending on ob- 
servation. In Fig.0we have plotted the average count rate for 
S a — 2.1 ± 0.1 versus S s hift- From Figs. [6] and 0a correla- 
tion between count rate and S'shift is observable: the two tracks 
with the lowest S'shift, tracks 1 and 8, have lower count rates 
for a given S a than the other tracks, except for the lowest three 
S a intervals of track 1 . Tracks 9 and 5 which have the highest 
Sshift show the highest count rates above S a ~ 2.5. 

3.3. The power spectra 

Power spectra with a Nyquist frequency of 2048 Hz were cal- 
culated using 16 s segments of high time resolution data and 
all energy channels, i.e., covering the 2-60 keV PCA band. 
Each of the 9 tracks separately was divided into (typically 8) 
S a intervals and the power spectra corresponding to the se- 
lected colour-colour points were averaged. A constant Poisson 
level determined by using the 1024 to 2048 Hz range (where no 
kHz QPOs were detected) was subtracted. Three representative 
power spectra from track 7 are presented in power x frequency 
vs. frequency representation in Fig. (HI The VLFN, BLN and a 
QPO are clearly visible. 

The power spectra were fitted with 2-4 Lorentzians. We 
chose to fit the VLFN with two zero-centered Lorentzians in- 
stead of a power law, because it has a "bumpy" structure, which 
can not be fitted well with a single power law (to obtain a good 
fit one has to add additional components to fit the bumps in 
any case) and is too broad a feature to be well fitted by a sin- 
gle Lorentzian. We will call these two components the low and 
high VLFN. This approach has the additional advantage that it 
avoids the common problem that the VLFN power law attempts 
to also fit power at high frequencies, which is probably unre- 
lated to the VLFN, and can bias the power law to ill represent 
the data at low frequencies. Two more Lorentzians fit the BLN 
and a QPO, when present. The Lorentzia ns are r eporte d in the 
i/ max , Q, r representation used by Belloni et al. (2002), where 
^max is the frequency at which the Lorentzian reaches its max- 
imum in a power x frequency vs. frequency power spectral rep- 
resentation, Q the quality factor (for which we use the standard 
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Fig. 5. The colour-colour diagram and hardness-intensity diagrams of track 1 . The colour-colour points have different symbols 
based on their location in the CD to simplify tracking the differences in the shape of the track. The count rates are for 5 PCUs. 
Typical error bars are shown for two points at the ends of the track in the different diagrams. 



definition of centroid frequency divided by full width at half 
maximum) and r the fractional rms amplitude. We will regard 
i/ max to be the characteristic frequency of a Lorentzian com- 
ponent. The rms powers were calculated in the 0-co Hz range, 
except for the two components fitting the VLFN where we used 
the 0.125-co Hz range because for these components the un- 
measured power below our lowest Fourier frequency provides 
a significant contribution (which is not the case for the other 
components). The rms amplitudes were calculated assuming 
an average background of 140 cnts s~ x (2-60 keV, 5 PCUs); 
the —20 cnts s _1 deviations from this average are negligible in 
view of the high source count rates. 

In all of the tracks we have found VLFN and BLN and in 
5 of the 9 tracks we de tected the 57-69 Hz QPO first reported 
bv lHoman etall |l?98). For all 5 a intervals we found VLFN; 
the BLN usually disappeared in the highest 5 a intervals and 
QPOs were only detected in the lowest 5 a intervals. The power 
spectral components and their dependence on 5 a and 5 s hift is 
discussed in more detail in Sect. 13.41 

3.4. Dependence of the power spectral properties 

on 5 a and S eh m 

3.4.1 . The very low frequency noise 

As described in Sect. 13.31 the VLFN was fitted with two 
zero-centered (Q = 0) Lorentzians. The power of both these 
Lorentzians was added to calculate the total rms amplitude of 
the VLFN. For all 9 tracks the rms of the total VLFN generally 
decreases with S a (see Fig.|9}, from -4.5-5.5% to -3-3.5%. 
An anti-correlation with S^ift is visible: tracks 5 and 9, which 
have a high S'shift, have a VLFN rms that is significantly lower 
than that of the other tracks. To quantify this we fitted a straight 
line to the VLFN amplitude - S a relation, and determined the 
value of the VLFN amplitude at 5 a = 2 from this fit. In Fig. 
EI| we show this amplitude as a function of S^ift- As 5 s hift 
increases from to —0.3 the VLFN rms at 5 a = 2 decreases 
from -4.8% to -3%. 



The absolute rms of the VLFN as a function of S a of all 
tracks combined is shown in Fig.[^ Up to S a = 2 the absolute 
rms decreases from —50 to 30 cnts s _1 , and for > 2 it 
stays approximately constant at —30 cnts s _1 . So for S a < 
2 the decreasing absolute rms is correlated to the decreasing 
count rate; however, for 5 a > 2 the count rate in some tracks 
increases, but this does not affect the absolute rms, which stays 
approximately constant. 

The characteristic frequency (i/ max ) of the low and high 
VLFN does not show a clear correlation with S a , but varies 
between 0.1-0.3 Hz for the low and between 0.5-1.5 Hz for 
the high VLFN (with a few exceptions up to 2 Hz when no 
BLN was fitted, and for low S a in track 9 where the BLN also 
shows an extremely high f max ). 

3.4.2. The band limited noise 

The BLN characteristic frequency, shown in Fig. [2] and rms, 
shown in Fig. ^] both decrease with S a . For most tracks the 
frequencies decrease from —10 to —4-5 Hz, except for track 9 
which has a BLN frequency of —30 Hz for its lowest two S a 
intervals (S a -1.5-1.7). The highest BLN rms is -4-6% and 
is found for low S a ; when S a increases the BLN rms decreases 
until it becomes undetectable around S a —2.5-3. Typical 95% 
confidence upper limits for the BLN for the S a intervals where 
it was not detected are —2%, with i/ maK = 5 Hz and Q = 0.5. 
For tracks 9 and 4 the rms increases again with 5 a above S a — 
3. 

The characteristics of the BLN do not show a clear corre- 
lation with Sshift, although the BLN in track 9 has a higher 
rms for a given S a than in the other tracks and also shows ex- 
tremely high frequencies for low S a (see above). The Q of the 
BLN is also not significantly correlated with S a . It varies be- 
tween and 1 , with the exception of the last interval of track 2 
(5 a = 3.13 ± 0.23, Q = 2.6tf 7 1 ) and S a interval 7 of track 
9 (5 a = 2.94 ± 0.05, Q = 47:*), where the BLN also shows a 
very high frequency of 12.1+Qg Hz; in these two cases (only) 
this component is sufficiently narrow to be called a QPO. 
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3.4.3. Quasi-periodic oscillations 

QPOs were detected with frequencies ranging from ^60-75 Hz 
and Q values of ^2-5 (detailed properties are shown in Table 
|2ji. They are only detected at low S a values, with the highest 
rms at the lowest S a , and in too small a range of S a to say 
anything about the dependence of their frequency and Q on 
S a . The fact that QPOs were detected in tracks 1 and 8 and not 
in tracks 5 and 9 could be related to the high S s hiii values of 
tracks 5 and 9, but this could also be due to the fact that these 
tracks do not reach the upper-left most part of the track, where 
the QPOs are found. The frequencies of the detected QPOs are 
all consistent wi t h bein g in the range of those previously found 
hv lHomanetal ] lll998h . 

We have searched for kHz QPOs up to 2048 Hz, but none 
were detected above the 3a confidence level. The strongest nar- 
row features (Q ~ 10) above 300 Hz were between 2 and 
3ct and resulted in 3a upper limits of up to 3% rms. A broad, 
3.1a, high frequency feature (2.6 ± 0.5 % rms, Q = 0.9±° 9 6 , 
v max = 858±T^ Hz) is found in track 7 (S a = 1.50 ± 0.08), 
but when the number of trials is taken into account this is not a 
significant detection. 
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Fig. 7. The count rate at 5 a = 2 for all 9 tracks plotted vs. 
<5 s hift- The two tracks with the lowest 5 s hift (tracks 1 and 8) 
show the lowest count rates. 
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Fig. 8. Representative power spectra of GX 13+1. Shown are power spectra of track 7 with an averaged S a of 1.50 ± 0.08 (left; 
x 2 /dof = 125/145), 2.01 ± 0.06 (middle; x 2 /dof = 135/148) and 2.73±0.12 (right; x 2 /dof = 138/151). All power spectral 
features become weaker with increasing S a . The solid line represents the fit and the dashed and dotted lines the Lorentzian 
components. The VLFN is fitted by the long dashed line (low VLFN) and the dotted line (high VLFN), the BLN by the dash- 
single dot line and the QPO by the dash-triple dot line. In the left power spectrum some excess power is present in the 200-2000 
Hz range, but no significant kHz QPO is detected. 

track 1, 12, 4, 6 and 9 track 8, 2, 7 and 5 




Fig. 9. VLFN fractional rms amplitude vs 5 a for all 9 tracks; the amplitude plotted is that of the combined low and high VLFN, 
integrated from 0.125 to oo Hz. For all tracks the amplitude of the VLFN generally decreases with S a , and tracks 5 (five-pointed 
stars, right frame) and 9 (six-pointed stars, left frame), which have the highest S s hift, have the lowest amplitude for a given S a . 
Symbols are the same as in Fig.|6] 
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Fig. 10. The VLFN amplitude at S a =2 for all 9 tracks (as deter- 



mined from a linear fit to a plot of the VLFN rms vs. S a ), as a 
function of 5 s hift. The amplitude decreases from ~4.8% to 
as Sghift increases from to ~0.3. 



Fig. 11. The absolute amplitude of the VLFN as a function of S a 



,y for all 9 tracks combined. Up to S a =2 the absolute rms ampli- 
tude decreases from ^50 to 30 cnts s~\ and for S a >2 it stays 
approximately constant at ^30 cnts s _1 . 
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2 3 4 1 2 3 



Fig. 12. The characteristic frequency of the BLN vs S a for all 9 tracks. The characteristic frequency generally decreases with 
S a . Track 9 shows an extremely high frequency for its lowest two S a intervals. Symbols are the same as in Fig.|6j in some S a 
intervals no BLN was detected. 

track 1, 12, 4, 6 and 9 track 8, 2, 7 and 5 




2 3 4 12 3 



a a 

Fig. 13. The fractional rms amplitude of the BLN vs S a for all 9 tracks. The amplitude of the BLN generally decreases with S a 
for all tracks, but in track 4 (diamonds, left frame) and 9 (six-pointed stars, left frame) it increases again for the highest two S a 
intervals. 3cr Upper limits were determined with Q = 0.5 and v max = 5 Hz, and are represented by the points with only a a 
negative error bar extending to the bottom of the plot. Symbols are the same as in Fig. [6] 
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4. Discussion 



4.3. GX 13+1 in the Z/atoll scheme 



4. 1. Extreme secular motion in CD and HIDs 



We have studied the secular motion both in the CD and in the 
HIDs of GX 13+1. The largest shifts are observed in the direc- 
tion of HC (-40%) and not in SC (-7%). In other atoll and Z 
sources which show secular motion, shifts of less than 3% in 



Di Salvo et alJ 


2003) 


Kuulkers et al. 


1994) 



are 



6% in HC and 8% in SC (GX 5-1 
observe d; in Cyg X-2 the shifts in both HC and SC are up to 
-10% JKuulkerset alJ ll996). This is not strongly dependent 
on the choice of bands. The shifts in the CD of GX 13+1 are 
remarkably large, particularly in HC, even when compared to 
Z sources, where the shifts are larger than in atoll sources. As 
is also seen in other sources (references above), the shift of the 
track shows a correlation with average count rate resulting in 
shifted tracks in the HIDs. We find that in GX 13+1 the am- 
plitude of the VLFN also shows a correlation with this shift; 
it decreases when the track is shifted towards higher HC and 
SC. Such a correlation is not observed in the Z source GX 5- 
1, but in the atoll source 4U 1636-53 the part of the track that 
shifted to higher SC and intensity also has a lower VLFN rms. 
The VLFN rms in Cyg X-2 shows the same dependence on 
the shift of the track, but the overall intensity level of a track 
decreases with this shift jKuulkers et alj fl9991 In GX 13+1, 
track 9, which has the largest shift also has unusual BLN prop- 
erties (Sect. l3~4~2l i: in 4U 1636-53 an extra, broad component 
(Q = 0.8) is detected next to the BLN in the high intensity part 
of the track at — 1.1% rms, but not in the low intensity part 
(with la upper limits of —0.6%). 



4.2. The flipping effect 



In the HIDs with the HC or SC on the vertical, and count rate 
on the horizontal axis the arc-shape found in the CD is re- 
versed, i.e., the source moves through the arc in the opposite 
sense in the HIDs as compared to the CD. This is different from 
what is usually observed in atoll sources JWiinands et alJl999t 
van Straaten et alJ l2000t iMuno et alJ 120021; iDi Salvo et alJ 



2003), where, on short timescales, the same basic pattern is 
traced out in both the CD and the HIDs and the sense of mo- 
tion is always the same. On longe r timescales secular shif ts in 
the HIDs are sometimes observed (van Straaten et al. l2003h . In 
most Z sources t he observed basic pattern in the CD and HIDs 
is also similar llSchulzet alJll989t I Kuulkers & van derKlisI 
Il996l IJonker et alJ l2002HHoman et alJ 120021) . although for 
some sources the count rate decreases along (a part of) the FB, 
causing this branch t o be in the opposite d irection in the HIDs 
compared to the CD ( Kuulk eTs et alj fl996). The cause for this 
unusual behaviour of GX 1 3+ 1 in the HIDs is that the count rate 
depends on S a in a different way than in most other sources (see 
Sect. 13. 21 . A possible explanation for this will be discussed in 
Sect.l4~4l 



The tracks that GX 13+1 traces out in the CD are, at first sight, 
very reminiscent of a classic island/banana in an atoll source, or 
a normal-branch/flaring-branch in a Z source. However, from 
the results presented in Sect. [5] we find that the power spectral 
properties of this source as a function of S a resemble neither 
those of a classic atoll, nor of a Z source. 

In a given track the count rate of GX 13+1 decreases along 
the track up to S a =2-3 after which (in some tracks) it increases 
again with S a , while in atoll sources on short timescales the 
count rate generally increases monotonically with S a along 
the whole trac k (see Se ct. 13.21 f or refe r ences on atoll source 
behaviour see Ivan Straaten et all 120021; IDi Salvo et alJ 120031: 
Hasinger & van derKlisll 989). This is also different from what 
is observed in Z sources, where the count rate can slightly 
decrease along the NB, but increases immediately after the 
NB/FB vertex so that the minimum in count rate is at the ver- 
tex (see Sect. 13.2} for references on Z source behaviour see 



Dieters & van der KlislEoOol IJonker et alj[ 2000: O' Neill et alJ 
20011 IJonker et alJl2002t iHoman et alJl2002i) . The VLFN of 
GX 13+1 is very strong (usually —4-5.5% rms for S a <2) 
compared to atoll sources in the island state or lower banana 
branch and Z sources in the NB (the VLFN frequency range 
is similar). The VLFN rms of GX 13+1 generally decreases 
with S a . This is the opposite of what is normally seen in both 
atoll and Z sources, where the VLFN increases with S a . along 
the whole track (see, however, Wiinands et al. 1999, on the 
atoll source 4U 1820-30). The BLN is quite weak (-2-5%) 
compared to that in atoll sources in the island state and lower 
banana and its rms decreases with S a as is usual in atoll and 
Z sources, but in tracks 4 and 9 it increases again in the high- 
est S a intervals. The characteristic frequency of the BLN also 
decreases with S'a, which, in either atoll sources in the island 
state and lower banana branch or Z sources, has so far only 
been observed in the atoll source 4U 1820-30, where a BLN- 
like c omponent that evolve s into or is replaced by a —7 Hz 
QPO ( Wiinan ds et alJl999t) is seen. In the Z source GX 349+2, 
a similar noise component is seen of which the centroid fre- 
quency decreases with S'a, but the v max of this component 
shows a more complicated dependence on S a llO'NeilletalJ 
120021) . Generally, in atoll sources the ^ max of the BLN in- 
creases with S a up to the upper banana, where it starts to de- 
crease again, whereas in Z sources it increases monotonically. 

We find the QPO properties of GX 13+1 to be unlike ei- 
ther Z or atoll sources. At th e lowe st S a intervals the 57-69 Hz 
al] 1 1998) 

is detected. In atoll sources 



QPO found by Homanj 
a low frequency QPO is usually detected in the is land state with 
a frequency of 10-15 Hz (see, however, IWiinands et alj fl999) 
increasing up to 40-47 Hz towards the upper banana, where it 
becomes undetectable (although a second low frequency QPO 
is som etimes also seen at sim ilar frequencies, but for higher 
S a , see Ivan Straaten et al.l2002l) . If the QPO would remain de- 
tectable, one would expect to find such a 57-69 Hz QPO at 
higher S a , in the middle banana branch. The HBO in Z sources 
has a frequency up to —55 Hz in the upper NB; if the QPO in 
GX13+1 is interpreted as this QPO one would also expect to 
find a -6 Hz NBO on this branch; in the Z source GX 349+2 
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no N BO is found, but a very broad noise co mponent is seen in- 
stead jKuulkers & van der Klisll998UO'Neill et all200lb . and 
of this GX 13+1 shows no evidence either. The fact that no kHz 
QPOs are detected anywhere along the branch (see Sect. 13.4.31 
is uncharacteristic for both atoll and Z sources. 

Clearly, GX 13+1 is a peculiar source in which the proper- 
ties correlate to S a in a very unusual way when compared both 
to atoll and Z sources. In the final section we explore a possible 
explanation for this unusual behaviour. 

4.4. A jet scenario 

If the behaviour of GX 13+1 is interpreted within the frame- 
work of atoll source behaviour, then clearly for high S a the 
source is in the upper banana. Based on the CD, at low S a one 
would expect the source to be in the lower left banana and is- 
land state. However, the fact that for low S a 57-69 Hz QPOs, 
and weak BLN with a characteristic frequency that decreases 
with S a are found, suggests that, contrary to what one would 
expect based on the CD, the "true" state of the source there 
is the middle part of the banana branch, because that is where 
other atoll sources can be expected to show this kind of power 
spectral behaviour. This is supported by the fact that so far no 
kHz QPOs have been found in this source, consistent with their 
absence in other atoll sources in the middle and upper parts of 
the banana branch. 

In GX 13+1 the count rate decreases as S a increases up to 
S a =2— 3, and only for higher S& increases, which is what in a 
regular atoll source one would expect along the whole track 
(see Sect. 13.21 . One interpretation for this unusual behaviour 
would be that there is something adding extra flux at low S a , 
which becomes weaker towards higher S a resulting in a count 
rate that first decreases, and then increases as the standard atoll 
upper banana behaviour takes over. From the two branched 
structure in the CDs one can conclude that this extra flux would 
have to be hard, because for low S a the source spectrum gets 
harder when the count rate increases. This scenario would also 
explain why the BLN is so weak and shows a decreasing f max 
along the track, and why no kHz QPOs are detected at low S a : 
these properties are normal for the middle and upper parts of 
the banana branch, which is, by hypothesis, the true state of 
the source. In most atoll sources the low frequency QPO disap- 
pears at the same S*, or a slightly lower value, as the kHz QPO 
(Ivan Straaten et alJl2Q02t f pi Salvo et af]l2003l) : the fact that in 
GX 13+1 we still see the low frequency QPO while the kHz 
QPO is absent is unusual, but in line with the relatively high 
frequency of the low frequency QPO. So the hypothesis we are 
exploring is that the hard extra flux creates an anomalous sec- 
ond branch in the CD reminiscent of a lower banana and island 
branch, while the source remains in the "middle" banana state. 

A possible source for this extra, hard flux could be a rela- 

tivistic jet. In general in atoll and Z sources radio emission, as- 

j II 

cribed to the presence of a relativistic jet (Fender 2002), tends 
to increase towards lower S a . If this more intense jet is associ- 
ated with the extra hard flux, this could cause a second branch 
in the CD, while the source is still in the upper banana state. 
We further hypothesize that the unusually strong VLFN would 



occur in the hard jet flux. Consequently, when the jet becomes 
weaker towards higher 5 a this causes the VLFN to decrease in 
amplitude. This is in accordance with the fact that up to S a =2 
the absolute rms amplitude of the VLFN decreases, from which 
it can be concluded that it is the flux, that also decreases in this 
part of the track, that shows relatively strong VLFN that dis- 
appears. For S a >2 the absolute rms amplitude of the VLFN 
stays approximately constant, while in some tracks the flux in- 
creases for S a >3. The extra flux in these tracks does lead to a 
significant increase of the absolute rms amplitude. 

The question of course is: why would a jet in GX 13+1 
cause this peculiar behaviour and not in other atoll sources? 
The radio emission of GX 13+1 is unusually strong for an atoll 
source JFender & Hendrvll2000h . so perhaps for a reason un- 
known the jet is unusually strong in this source. A more at- 
tractive possibility is that it is not the jet that is unusual, but 
the system's orientation. If the jet is pointing directly towards 
us Doppler boostin g could si gnificantly increase its radio flux 
(e.g. Fender & Hendrv 2000); a strongly relativistic jet (which 
has already been d irectly imaged for the Z source Sco X-l, 
Fomalont et al. 2001a b) could easily Doppler boost the radio 
emission by more than an order of magnitude, putting the radio 
emission of GX 13+1 in accordanc e with measurements an d 
upper limits of other atoll sources JFender & Hendrvl fcoOO). 
The unusual X-ray phenomena could then for example be the 
consequence of a better view of the X-ray emitting regions as- 
sociated with the base of the jet, offered by the unusual ori- 
entation of the disk, possibly assisted by Doppler boosting as 
well. A similar possibility of Doppler-boosted jet X-ray emis- 
sion causing hard flarin g ("microblaz ar") was recently explored 
for Cvg X-1 hv lRomero et all J2002h . 

The large shifts of the track in the CD of GX 13+1 com- 
pared to other atoll and Z sources might be related to this better 
view of the X-ray emitting base of the jet region as well. This 
would require this region to contribute different baseline levels 
of hard flux at different epochs in order to produce the secular 
motion (which, consistent with this, takes place mostly in the 
HC), VLFN whose amplitude anticorrelates, as well as extra 
flux at S a > 3 which correlates, with this baseline flux level. 
This could for example be related to secular geometry changes 
(e.g. precession) modulating our view of the base of the jet, or 
to variations in the jet activity itself. 

One can also attempt to explain the peculiar behaviour of 
GX 13+1 with a similar relativistic jet scenario assuming the 
source is a Z source. The QPO observed in GX 13+1 would 
then be interpreted as the HBO, occurring at the top of the NB, 
leading to the conclusion that the CD shows both a NB and 
a FB. The strong VLFN that decreases along the whole track 
might be related to a jet similarly to the atoll scenario. This 
jet could be related to the high radio over X-ray luminosity 
ratio compared to other Z sources. However, this model leaves 
unexplained the decreasing frequency of the BLN with S a and 
the absence of an NBO. 

So, when a choice is made between an atoll or a Z clas- 
sification for GX 13+1 (in either case modified by an unusual 
jet geometry) the overall phenomenology favours the former 
option. An alternative would be that GX 13+1 combines some 
properties of atoll and Z sources (as noted in Sect.Q GX 13+1 
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has been described as the atoll source closest to the Z sources). 
This possibility wil l be further explored in a forthcoming paper 
faoman et all2003fc . 

It has been sug gested that narrow absorption features, as 
found in GX 13+1 (Ue da et al.ll200 1). are related to a high in- 
clination of the source, contrary to the low inclination predicted 
by the scenario discussed above. However, it is also possible 
that the absorption features are not related to the inclination o f 
the source, but, as previously suggested bv lUeda et alJ {l998), 
to a highly ionized plasma associated with their relativistic jets, 
causing absorption from the flux from the hot inner accretion 
region. 

Several predictions can be made based on the atoll 
source+jet hypothesis. Since the unusually strong VLFN oc- 
curs in the jet flux, one expects the VLFN, at low S a , to be the 
strongest in the hard flux. Also, as the jet is the main source 
of radio emission in LMXBs, we expect the radio emission to 
be strongest in the upper part of the second branch in the CD 
(at low Sg), decreasing as the source moves towards the vertex. 
This can be tested with simultaneous radio and X-ray obser- 
vations. One would have to account for a d elay factor in the 
radio emission of perhaps up to half an hour ein-Wolt et al. 

2002) associated with the expansion o f the synchro t ron clo uds. 
Such observations have been done bv lGarcia et alJ dl988t) . but 
with a small time range (~5 hours of overlapping X-ray and 
radio observations) and no CD analysis was performed, and by 
iBerendsen et alJ ll200Cl) . also with a short period of overlapping 
X-ray and radio observations. Recent simultaneous radi o and 
X-ray observations also exist (to be reported in iHoman et alJ 

2003) , and preliminary results of this work that became avail- 
able after the current analysis and interpretation were com- 
pleted, appear to confirm our suggestion. 
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Table 1. The RXTE observation IDs of the observations that are included in each track, the HC and SC of their vertices, and the 
distance, S^ift, over which each track is shifted relative to track 1. Average flux at the vertex (S a = 2), start time of the first and 
end time of the last included observations and the total included observation time are also indicated. Note that tracks 3, lOa-c and 
1 la-b do not show a vertex. All dates are in the year 1998; typical errors in HC are ~ 0.01, in SC^ 0.03, in S s hat ~ 0.01 and in 
flux at = 2 - 43. 
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Table 2. All detected QPOs and their properties. Typical 95% 
confidence upper limits (^ max = 68 Hz, Q = 3) for tracks in 
which the QPO was not detected are up to ~2% rms. 



